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Massive neutron stars may harbor deconfined quark matter in their 
cores. I review some recent work on the microphysics and the phenomenol- 
ogy of compact stars with cores made of quark matter. This includes the 
equilibrium and stability of non-rotating and rapidly rotating stars, gravita- 
tional radiation from deformations in their quark cores, neutrino radiation 
and dichotomy of fast and slow cooling, and pulsar radio-timing anomalies. 

PACS numbers: 97.60.Jd, 26.60.Kp. 95.30.Sf 

1. Introduction 

Because of the quark substructure of nucleons predicted by quantum 
chrome-dynamics (QCD), nuclear matter will undergo a phase transition to 
quark matter if squeezed to sufficiently high densities. In the quark matter 
phase the "liberated" quarks occupy continuum states which, in the low- 
temperature and high-density regime, arrange themselves in a Fermi sphere. 
The Fermi-sphere determines the form of low-lying excitation spectrum of 
quark matter, which resembles that of less exotic low-temperature systems 
found in condensed matter (e.g. electron gas or ultracold atomic vapor) 
or hadronic physics (e.g. nuclear or neutron matter). In analogy to these 
system the attractive interaction between quarks, mediated by the gluon 
exchange (which is responsible for the bound state spectrum of QCD, e.g., 
the nucleon and mesons) leads to quark superconductivity and superfluidity 
(color superconductivity) via the Bardeen-Cooper-Schrieffer mechanism pQ. 
Furthermore, under stellar conditions the pairing between the two light 
flavors of quarks occurs at finite isospin chemical potential, i.e., when the 
Fermi surfaces of up and down quarks are shifted apart by an amount which 
is of the same order of magnitude as the gap in the quasiparticle spectrum. 
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Under these conditions the pairing between the fermions persists, but the 
actual pairing pattern may be significantly different from that of the BCS 
and is likely to involve breaking of the spatial symmetries by the condensate 
order parameter. 

During the last decade there has been a substantial progress in under- 
standing of pairing in two-component asymmetric superconductors. Firstly, 
new developments on (variations of) the Larkin-Ovchinnikov-Fulde-Ferrell 
(LOFF) phase revealed novel lattice structures of the order parameter [2]. 
Secondly, it has been suggested that such systems may actually phase sepa- 
rate into normal and superconducting domains [3j or the pairing may require 
changes in the shapes of the Fermi surfaces [4|. Controlled experiments on 
two-component cold atomic vapors with mismatched Fermi surfaces show 
that the phase separation scenario is at work. Whether this is the case in 
the related systems such as the isospin asymmetric nuclear matter or decon- 
fined quark matter is still unclear. Therefore, a major goal of astrophysics 
of dense matter is to find and to quantify the manifestations of dense phases 
in the observable properties of compact stars. Here we review some recent 
progress towards this goal. 

2. Color superconducting phases with broken spatial symmetries 

The quark-quark interaction is strongest for pairing between up and 
down quarks which are antisymmetric in color pQ . At intermediate densities 
quark matter is composed of up and down quarks, while strange quarks can 
appear in substantial amounts at higher densities. The leading candidate 
phase at these densities is the so-called 2SC phase, which is characterized by 
the order parameter A oc {ip T (x)C^T2X2'4 , {x)} , where T2 is the Pauli matrix 
in the isospin state, A2 is the Gell-Mann matrix in the color space, C = 
Z7 2 7° is the matrix of charge conjugation. The isospin chemical potential 
is determined by the /3-equilibrium condition — \x u = \i e ~ 100 MeV 
among d and u quarks and electrons, where /x,; with i = d, u, e are the 
chemical potentials. The Fermi spheres of up and down quarks are shifted 
apart by this amount. Thus, the cross-flavor pairing must overcome the 
disruptive effect of mismatched Fermi surfaces. Furthermore, if the physical 
strange quark mass is close to its current mass ~ 100 MeV, electrons may 
be gradually replaced by strange quarks, which again will disfavor the cross- 
flavor pairing. 

A superconducting phase with asymmetric cross-species pairing needs 
to optimize the overlap between the Fermi surfaces (which is perfect in the 
symmetric BCS state) to attain the maximum possible condensation energy. 
There are a number of ways to achieve this: (a) condensate that carries non- 
zero momentum with respect to some fixed frame. The penalty in the energy 
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budget for the (always positive) kinetic energy of condensate motion is com- 
pensated by the gain in the (negative) condensation energy. In three-flavor 
quark matter the face-centered cubic (foe) lattices were identified as having 
particularly low free energy [2] in the regime where the Ginzburg-Landau 
(GL) theory is applicable. The spatial form of the condensate beyond the 
GL regime is not known, therefore current treatments beyond this regime 
require some simplifying assumptions [5], [6] . (b) A more symmetric super- 
conducting phase exploits deformations of Fermi spheres as the mechanism 
of restoring the coherence needed for pairing at the cost of kinetic energy 
loss caused by these deformations. The deformed Fermi surface phase re- 
quires minimal breaking of spatial symmetry and is more symmetric than 
the lattice phases above [4J. 

In this review we will be concerned mainly with the 2SC and the crys- 
talline color superconducting (CCS) phases. We expect that the phenomenol- 
ogy of related phases with cross-flavor pairing is essentially the same. 

3. Equilibrium and stability 

The central question, which we address in this section, is whether the 
equation of state of matter at high densities admits stable configurations 
of self-gravitating objects in General Relativity featuring deconfined quark 
matter. The deconfinement phase transition from baryonic to quark matter 
leads to a softening of the equation of state which could lead to an insta- 
bility towards a collapse into a black hole. In ref. [7] the quark phase was 
studied in the Nambu-Jona-Lasinio (NJL) model, which is a low-energy non- 
perturbative approximation to QCD, that is anchored in the low-energy phe- 
nomenology of the hadronic spectrum. While dynamical symmetry break- 
ing, by which quarks acquire mass, is incorporated in this model, it lacks 
confinement. Physically, the true nuclear equation state must go over to 
some sort of quark equation of state at some density if the deconfinement 
has to take place in nature. If pressure vs chemical potential curves for 
chosen equations of state of nuclear and quark matter do not cross, then 
the models are incompatible in the sense that they can not describe the 
desired transition between the nuclear and quark matter. The low-density 
equation of state of nuclear matter and the high-density equation of state 
of CCS matter were matched in ref. [7J at an interface via the Maxwell 
construction. The phase with largest pressure is the one that is realized 
at a given chemical potential. Thus, at the deconfinement phase transition 
there is a jump in the density at constant pressure as illustrated in Fig. [H 
left panel. The low-density nuclear equations of state are based on the 
Dirac-Bruckner-Hartree-Fock approach and are the most hard equations of 
state in our collection [9]. It should be noted that the matching with 
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Fig. 1. Left panel. Number density versus pressure for three models. For the 
models A and Al the nuclear (low density) equation of state is the same; for the 
models A and B the quark (high density) equation of state is the same. Right 
panel. Dependence of the total stellar mass and the mass of the quark core in units 
of solar mass Mq on the central density for non-rotating configurations. The lower 
set of curves represents the masses of the CCS quark cores, the upper set - the 
total masses of the configurations. The maximal masses are marked with boxes. 



softer equations of state can be enforced by varying the normalization of 
the pressure ("effective bag constant") of the quark matter, in which case 
the density of deconfinement is a free parameter. 

The spherically symmetric solutions of Einstein's equations for self- 
gravitating fluids are given by the well-known Tolman-Oppenheimer-Volkoff 
equations [8 J . A generic feature of these solutions is the existence of a max- 
imum mass for any equation of state; as the central density is increased 
beyond the value corresponding to the maximum mass, the stars become 
unstable towards collapse to a black hole. A criterion for the stability of 
a sequence of configurations is the requirement that the derivative dM/dp c 
should be positive (the mass should be an increasing function of the central 
density) . 

For configurations constructed from a purely nuclear equation of state 
the stable sequence extends up to a maximum mass of the order 2 Mq (Fig. 
[1]), this large value being a consequence of hardness of the equation of state. 
The hybrid configurations branch off from the nuclear configurations when 
the central density reaches that of the deconfinement phase transition. It 
is seen that a stable branch of hybrid stars emerges in the range of central 
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densities 1.3 < p c < 2.5 x 10 15 g cm" 3 . The masses of the CCS quark 
cores cover the range < M core /MQ < 0.75 — 0.88 for central densities 
1.3 x 10 15 < p c < 2 x 10 15 . Thus, the quark core mass ranges from one third 
to about the half of the total stellar mass. 

Millisecond neutron stars can rotate at frequencies which are close to 
the limiting orbital Keplerian frequency at which mass shedding from the 
equatorial plane starts. The Keplerian frequency sets an upper limit on the 
rotation frequency, since other (less certain) mechanisms, such as secular 
instabilities, could impose lower limits on the rotation frequency. The mass 
versus central density dependence of compact stars rotating at the Keplerian 
frequency is similar to that for non-rotating stars with the scales for mass 
shifted to larger values [7j. The increase in the maximum mass for stable 
hybrid configurations (in solar mass units) is 2.052 — ► 2.462 for the model 
A, 2.017 — ► 2.428 and 2.4174 for the model Al (there are two maxima) and 
1.981 -> 2.35 for the model B [7]. 

Thus, a new branch of stable hybrid configurations with CCS matter 
cores emerges within a broad range of central densities. The quark equation 
of state depends only on the NJL model parameters (i.e. does not contain an 
additional "bag constant" ) . A quite general conclusion of our analysis is that 
the nuclear equation of state needs to be hard to enable thermodynamics! 
equilibrium with the NJL model quark matter. 

4. Gravitational radiation 

Continuous gravitational waves emitted by non-axisymmetric rotating 
compact stars are expected to be in the bandwidth of current gravitational 
wave interferometric detectors. Upper limits on the strain of gravitational 
waves from a selection of radio pulsars were set by the LIGO collabora- 
tion |10| . Gravity waves arise from time-dependent quadrupole deformations 
of masses. Therefore, rotating isolated neutron stars will emit gravitational 
radiation if their mass distribution is non-axisymmetric with respect to their 
rotation axis. The axial symmetry of the star's core can be broken by the 
solid deformations in CCS matter, whose shear modulus ^ s hear was com- 
puted in ref. [TT]. Initially, the gravitational wave emission from a quark 
star made of uniform-density incompressible CCS matter was estimated by 
Lin [12]. Haskell et al. [12J estimated core deformations and associated 
with them constraints on the QCD parameters for sequences of 1.4 Mq and 
R = 10 km stars composed of CCS incompressible quark cores and hadronic 
shells obeying n = 1 polytropic equation of state. Gravitational radiation of 
models based on realistic equations of state were considered in ref. [13] . One 
key difference to the previous models is that realistic hybrid configurations 
have masses that are close to the maximum sustainable mass M max ~ 2Mq. 
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Fig. 2. Left panel. Dependence of density on the internal radius (model A). The 
dashed lines correspond to purely nuclear, the solid lines - to hybrid configurations. 
Density jumps at the phase transition to quark matter. Right panel. The strain of 
gravitational wave emission of the Crab pulsar. Each triple of curves corresponds 
to the gap value (from top to bottom) A = 50, 40, 30, 20, and 10 MeV and the 
breaking strain <r max = 10~ 4 . The horizontal line shows the current experimental 
upper limit. 



Furthermore, the quadrupole moment of the quark core was computed from 
realistic density profile, which differs substantially from constant density 
profile of an incompressible fluid, see Fig. [2j The characteristic strain am- 
plitude of gravitational waves emitted by a triaxial star rotating about its 
principal axis is ho = 16ir 2 Gc~ 4 eI zz i' 2 r^ 1 , where v is the star's rotation 
frequency, r is the distance to the observer, e = (I xx — I yy )/I zz is the equa- 
torial ellipticity, Iij is the tensor of moment of inertia, G is the gravitational 
constant, c is the speed of light. The elastic deformations are assumed to 
be small perturbation on the background equilibrium of the star. Fig. 
displays the strain of gravitational wave emission as a function of the core 
mass M core and the current upper limit for the Crab pulsar ho < 3.4 x 10 -25 , 
which is rotating at the frequency v = 29.6 Hz at the distance 2 kpc [10j . 
The value of the breaking strain of CCS core is <7 max = 10~ 4 . The depen- 
dence of the stain amplitude on the microscopic parameters follows from 
the chain ho oc Q ma , x oc ^ s hear oc <r max A 2 , whereby Q ma , x is the maximal 
quadrupole moment of the CCS core. The strain of gravitational wave 
emission is close to the upper limit for A = 50 MeV. If, however, the gaps 
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are small more "optimistic" values of 10 3 < <7 max < 10 2 will be needed to 
generate sizeable strain amplitude. The Crab puslar's current limit implies 
<r max A 2 ~ 0.25 MeV 2 , assuming maximally strained matter. The evolution- 
ary avenues that may lead to such maximal deformations are not known. 
Pulsars may simply preserve their initial deformations as they cool down 
and solidify in the CCS state. We conclude that massive compact stars with 
CCS quark cores are strong candidate sources of gravitational wave emis- 
sion. Improved upper limits on the strain amplitude, can narrow down the 
admissible range of parameters of the CCS matter in the future. 

5. Neutrino emission and cooling 

Neutrino emissivities control the cooling rate of a neutron star during the 
first 10 4 — 10 5 yr of their evolution. At later times the photon emission from 
the surface and the heating in the interior are the main factors. Depending 
on the dominant neutrino emission process in the neutrino emission era 
t < 10 5 yr the cooling could be slow (standard) or fast (nonstandard). 

The slow cooling scenario is based on neutrino cooling via the modi- 
fied Urea and bremsstrahlung processes. The fast cooling requires uncon- 
ventional processes such as decays in pion/kaon condensates, direct Urea 
process on nucleons, hyperons, or quarks. Phase-space arguments show 
that fast cooling neutrino processes have typically temperature dependences 
~ T 6 , whereas those leading to slow cooling require additional phase space 
for the spectator particle and their emissivities scale as T 8 . 

An inspection of the observational data on neutron star surface tem- 
peratures, which is commonly presented on a plot of photon-luminosity (or 
surface temperature) vs age, shows that the data cannot be described by a 
single cooling track; a regulator is needed that will cool some stars faster 
than the others [9]. It is reasonable to assume that the heavier stars 
cool via some fast mechanism, while the lighter stars cool slowly via the 
modified processes. Indeed, the fast cooling agents operate above a certain 
density threshold. Therefore the dichotomy of observed high and low sur- 
face temperature could be an evidence of a fast process operating in massive 
stars. Thus, we are led to examine the potential fast cooling processes in 
the quark cores of compact stars. The key process is the Urea processes 
(quark /3-decay) d — > u + e + v e and u + e — > d + u e , which is permitted in 
interacting quark matter for any asymmetry between u and d quarks. The 
effect of asymmetric pairing is well demonstrated on the example of 2SC 
phase |14[ [To] . The result for the emissivity depends on whether the pa- 
rameter z = A/ 5/j, is larger or smaller than unity, where A is the gap in the 
2SC phase and Sfi = ^d — fJ- u is the shift in the chemical potentials of the up 
and down quarks. For z > 1 all the particle modes are "gapped" , therefore, 
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as the temperature is lowered, the emissivity is suppressed (for asymptoti- 
cally low temperatures exponentially). When z < 1 there are gap less modes 
in the quasiparticle spectrum; this implies that the neutrino production is 
not affected by color superconductivity. As a result, the superconducting 
quark matter cools at a rate comparable to the unpaired matter. More 
generally, we may conclude that the existence of ungaped segments on the 
Fermi surfaces of quarks will lead to an enhanced cooling of these stars. A 
density dependent z parameter may resolve the dichotomy of the fast and 
slow cooling by providing a smooth transition between these extremes. 

6. Radio-timing and rotational anomalies 

Pulsars are nearly prefect clocks, whose pulsed emissions, with a typi- 
cal periodicity of seconds or less, are locked to the rotation period of the 
star. Their periods increase gradually over time, due to a secular loss of 
rotational energy. Some pulsars show deviations from this regularity. The 
pulsar timing anomalies divide roughly into three types, (i) Glitches. These 
are distinguished by abrupt increases in the rotation and spin-down rates 
of pulsars by amounts AQ/n ~ 1CT 6 - 10~ 8 and AU/U ~ 1CT 3 . After a 
glitch, Af2/f2 and AQ/fl, slowly relax toward their pre-glitch values, on a 
time scale of order weeks to years, in some cases with permanent hysteresis 
effects. Such behavior is attributed to a component within the star that is 
only weakly coupled to the rigidly rotating normal component responsible 
for the emission of pulsed radiation, (ii) Timing Noise. These represent 
irregular, stochastic deviations in the spin and spin-down rates that are 
superimposed on the near-perfect periodic rotation of the star. Whether 
the superfluids are involved in the generation of timing noise remains un- 
clear, (hi) Long-Term Periodic Variabilities. Observed in the timing of few 
pulsars, most notably PSR B1828-11, these deviations strongly constrain 
theories of superfluid friction, if their periodicities are interpreted in terms 
of precession. (Note however that the Tkachenko modes of the vortex lattice 
are a viable alternative to the precession interpretation [17] ■) The impor- 
tance of the inferred precession mode stems from the fact that it involves 
non-axisymmetric perturbations of the rotational state, removing the de- 
generacy with respect to £ (£' ~ 0) that is inherent in the interpretation of 
post-glitch dynamics. 

The macroscopic physics of neutron-star rotation and its anomalies ob- 
served in the timing of pulsars can be described within the hydrodynamic 
theory of superfluids. At the local hydro dynamical scale, the rate of angu- 
lar momentum transfer between the superfluid and normal components is 
determined by the equation of motion of a vortex line in a neutral superfluid 



u s (v s -v L ) x v + ((v L - v N ) + C'(v L - V N ) X V = , 



(1) 
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where v$ and vjy are the superfluid and normal fluid velocities, vl is the 
velocity of the vortex, v is a unit vector along the vortex line, u>s is the 
unit of circulation, and £, £' are (dimensionless) friction coefficients, also 
known as the drag-to-lift ratios. These coefficients encode the essential 
information on the microscopic processes of interaction of vortices with the 
ambient unpaired fluid. 

The potential role of quark matter in the rotational dynamics of com- 
pact stars is not understood yet. Some of the candidate phases, such as 
the color-flavor-locked (CFL) phase, rotate by creating vortices [18] . which 
interact with the phonon gas within this phase [19J. However, phonons 
do not interact electromagnetically with the rest of the star, therefore it 
is unclear on which timescales the CFL quark core will couple to the ob- 
servable crust. Such a coupling is manifest for the 2SC and related phases, 
where the magnetic field (partially) penetrates into the superfluid region. 
While the vortices in these phases (if any) are not rotational in nature, their 
electromagnetic coupling to the vortices in the hadronic core and to elec- 
trons, which are shared with less denser phases of the star, can impact the 
rotational dynamics of pulsars. 

I would like to thank J. W. Clark, N. Ippolito, P. Jaikumar, B. Knippel, 
D. H. Rischke, C. D. Roberts, and F. Weber for collaboration on the topics 
discussed in this review. This work was in part supported by the Deutsche 
Forschungsgemeinschaft (Grant SE 1836/1-1) 
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